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Spin dynamics in pressure-induced magnetically-ordered phases in PHCC 
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We present inelastic neutron scattering experiments on the S = 1/2 frustrated gapped quantum 
magnet piperazinium hexachlorodicuprate under applied hydrostatic pressure. These results show 
that at 9 kbar the magnetic triplet excitations in the system are gapless, contrary to what was pre¬ 
viously reported. We show that the changes in the excitation spectrum can be primarily attributed 
to the change in a single exchange pathway. 
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I. INTRODUCTION 

Magnetic insulators offer a broad range of structural 
motifs, dimensionalities and well-defined short-ranged in¬ 
teractions. They are amenable to numerical modelling 
and are useful prototypes in the study of quantum phase 
transitions^ (QPTs). An important and well-studied case 
of a QPT is the so-called Bose-Einstein condensation 
(BEG) of magnons in gapped quantum magnets in ap¬ 
plied magnetic fieldsThese are soft mode transitions, 
where at the QGP excitations have a parabolic disper¬ 
sion, so the dynamical critical exponent z = 2. A quali¬ 
tatively different and much rarer soft mode transition can 
sometimes be induced in gapped spin systems through a 
continuous change of exchange constants. This, in turn, 
may in certain cases be achieved by the application of 
hydrostatic pressure. In these transitions the spectrum 
is expected to be linear at the quantum critical point 
(QGP), and hence z = 1. Until recently only one good 
experimental realization of a pressure-induced QPT had 
been found, namely that in the three-dimensional dimer 
system TlGuGls^^^— Eurther work has lead to fascinat¬ 
ing insights, in particular to the observation of a massive 
amplitude mode^i^ the magnetic analog of the Higgs bo- 
son4 

Quantum magnets built from organic molecules can 
be very susceptible to perturbation by external pressure 
due to their ‘soft’ molecular frameworksThis paper 
is concerned with the S = 1/2 quasi-two-dimensional 
gapped quantum antiferromagnet piperazinium hex¬ 
achlorodicuprate [(G 4 Hi 2 N 2 )Gu 2 Cl 6 , hereafter PHGG]. 
PHGG crystallizes in the triclinic space group PI with 
lattice parametersiS a = 7.984(4) A, 5 = 7.054(4) A, 
c = 6.104(3) A, a = 111.23(8)°, /3 = 99.95(9)°, 7 = 
81.26(7)°. The spin-1/2 Gu^+ ions are connected by 
a complex layered network (see Eig. El). The magnetic 
interactions in PHGG have been studied using inelastic 
neutron scatter ingi^ (INS) and were found to be highly 
frustrated and more complicated than a simple dimer 
model. The ground state is a spin singlet separated 
by a gap A = 0.98(6) meV from an S' = 1 triplet. 
More recent studies have demonstrated that the gap can 


be reduced by applied hydrostatic pressure, much as in 
TlGuGls. At p = 9 kbar the gap was found to decrease 
to A = 0.55 meV, and extrapolates to zero at ^ 20 kbar, 
hinting at a possible QPT at that point In contradic¬ 
tion with this result, more recent muon-spin relaxation 
(/i+SR) experiments have discovered that the destruc¬ 
tion of the spin-singlet state and the onset of magnetic 
long-range order occur at a much lower pc ~ 4.3 kbar^i^ 
A magnetically ordered phase of a Heisenberg spin sys¬ 
tem must have gapless spin waves, yet Ref. [ll| reported 
a spectral gap persisting well beyond the transition pres¬ 
sure observed in muon experiments. 

The purpose of the present inelastic neutron study is 
to resolve the apparent controversy. Our new, higher 
resolution experiments reveal that the spin excitation 
spectrum in PHGG is gapless in the pressure-induced or¬ 
dered state, already at 9 kbar. Moreover, we show that 
pressure-induced changes in the excitation spectrum can 
be primarily attributed to a variation in a single exchange 
pathway. The /i+SR experimentsi^ found two distinct 
magnetically ordered phases: an incommensurate phase 
above a critical pressure Pc ~ 4.3 kbar and a commensu¬ 
rate phase above pi ~ 13.4 kbar. Within the resolution 
of our experiment, we find no significant qualitative dif¬ 
ferences in the magnetic excitation spectra in the two 
regions. 

II. EXPERIMENTAL DETAILS 

High-quality deuterated single crystal samples of 
PHGG were grown by the thermal gradient method^i^ 
The sample was a single crystal of mass 150 mg. Before 
this experiment, the sample was successfully used in a dif¬ 
ferent neutron scattering study. The crystal was placed 
in a teflon tube and mounted with its 6 -direction vertical 
inside of a NiGrAl (“Russian alloy”) pressure cell. The 
top and bottom edges of the pressure cell (away from the 
sample) were covered in cadmium to reduce the back¬ 
ground from the cell. Eluorinert was used as pressure- 
transmitting medium. The pressure was applied in a 
piston press, measured mechanically at ambient temper¬ 
ature and then extrapolated to low temperature using 
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FIG. 1. (color online). False-color slices through the 4-dimensional inelastic neutron data sets measured in PHCC at T = 1.5 K, 
showing the dispersion of magnetic excitations in the (h, 0, 1) plane. The integration along the k and I or h directions is ±0.1 r.l.u. 
Each column corresponds to a different applied pressure. Top and middle: background-subtracted experimental data. Bottom: 
simulation, based on the model cross section and the fitted parameter values, and convoluted with the resolution, as described 
in the text. In all cases, the solid lines are the model magnon dispersion relation, plotted using the fitted parameter values. 


the documented behavior of the pressure cell. The pres¬ 
sure cell was mounted inside of a ^He cryostat on the 
time-of-flight cold neutron multi chopper spectrometer— 
(CNCS) at the Spallation Neutron Source (SNS), Oak 
Ridge National Laboratory. Data were collected with in¬ 
cident energy of E[ = 4.2 meV at T = 1.5 K. A converg¬ 
ing guide was used to focus the neutron beam vertically 
onto the sample. The spectra were recorded by mak¬ 


ing 180° rotations with 1° step size. Corrections for the 
energy-dependent transmission of the pressure cell were 
performed using a reference measurement of incoherent 
elastic scattering from a plastic test sample. Recorded 
events were projected onto the sample’s reciprocal space 
coordinate system and binned into two-dimensional cuts 
with 150, 150, and 80 bins for /i, /, and Tiuo^ respectively, 
covering the range (—3.7 r.l.u., —2.8 r.l.u., —0.5 meV) to 
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(3.7 r.l.u.,2.8 r.l.u.,3.5 meV). The data were integrated 
along k by ±0.1 r.l.u. and exported for further analysis 
using the MANTIEU^ program. 


at 9 kbar, with —0.4 < k < —0.2 is shown in Fig. [2j 
Note that at h = 1.5 the excitation energy has increased 
to ^ 0.7 meV, as compared to ^ 0 for /c = 0 in the 
top-center panel of Fig. [TJ 


III. RESULTS 

Our first result is that the lattice parameters a, c, and 
P do not change abruptly between 0 and 18 kbar indi¬ 
cating the absence of a structural transition. From the 
position of the Bragg peaks in the scattering plane we 
found: a = 7.8(8) A, c = 6.0(6) A, /3 = 99.9(9)° for 9 
kbar and a = 7.8(8) A, c = 6.0(6) A, ^ = 100.2(2)° for 
18 kbar (both at 1.5 K). 

Selected inelastic spectra taken at 0, 9, and 18 kbar 
are shown in Fig. [H These are two-dimensional slices 
through the four-dimensional data set, integrated along 
the k and h (or 1) directions in the range ±0.1 r.l.u. A 
smooth non-magnetic background was subtracted. For 
each cut, at each pressure, it was estimated by a second- 
order polynomial fitted to the intensity in areas where 
magnetic excitations are clearly absent. Even by eye it is 
immediately apparent that the excitations at 9 kbar are 
gapless with zero energy at the antiferromagnetic zone- 
center (1.5,0,—0.5). This point coincides with the mag¬ 
netic propagation vector of the ordered state that can 
be induced in PHCC at ambient pressure by an external 
magnetic field through a BEC-like transition^i^ 


IV. ANALYSIS 



FIG. 3. (color online). Structure of PHCC looking onto the 
crystallographic (a, c) plane. Some bonds have been omitted 
for clarity. The relevant Cu-Cu super-exchange pathways are 
indexed for the discussion in the text. 
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FIG. 2. (color online). False-color slice through the 4- 
dimensional inelastic neutron data sets measured in PHCC 
at p = 9 kbar, T = 1.5 K, showing the dispersion of magnetic 
excitations along the {h, —0.3, —0.5) direction. The solid line 
is a guide for the eye. 


Our setup, with the crystallographic 6-axis mounted 
vertically, and with a vertically-focusing guide, is ill- 
suited for the study of the dispersion of excitations along 
the 6* direction. Correspondingly, all the data shown in 
Fig. [Dare narrow slices around the (h,0,/) plane. All 
the analysis given in the next section is applied only to 
the data taken in that plane. This is important to keep 
in mind, since the limited data available does indicate a 
small but significant dispersion along 6*. One such cut 


Our analysis closely follows that used by Stone et al. 
in Ref. [H To within experimental resolution, the scat¬ 
tering can in all cases be attributed to a single sharp 
magnon branch. Correspondingly, the fitting model for 
the dynamic spin structure factor was written as: 

5'(Q, uj) = S{Q)5{nuj -hujQ). (1) 

Note that we are not making any distinction between dy¬ 
namic structure factors of different polarization, assum¬ 
ing spin correlations to be isotropic. Assuming PHCC 
is a Heisenberg spin system, this approach is fully justi¬ 
fied at ambient pressure, where the ground state is dis¬ 
ordered. However, in the magnetically ordered state the 
rotational symmetry is spontaneously broken. As a re¬ 
sult, correlations between different spin components, as 
well as the corresponding contributions to the dynamic 
structure factor, will no longer be equal. Nevertheless, 
considering the resolution and noise level of our data, 
we applied the isotropic model [Eq. ([1])] as an empirical 
function at all experimental pressures. 

In our model, also in line with Ref. IS the equal-time 
structure factor 5'(Q) was written as: 

5(Q) oc Yi;d(cos(Q-d)-l). (2) 

TlLdQ 

In this expression the sum runs over all pairs of spins 
with relevant exchange interactions, separated by a bond 
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d. Once again, as long as the system is isotropic, this 
expression is well justified. It directly follows from the 
single-mode approximation [Eq. ([T])], and the Hohenberg- 
Brinckman sum rule^^ for the first moment of the dy¬ 
namic structure factor. As explained in detail for ex¬ 
ample in Ref. [S the coefficients in this case are 
the expectation values of the Heisenberg exchange en¬ 
ergy on the corresponding bonds between spins: Ed = 
Jd(So • Sd). In the magnetically ordered (anisotropic) 
phase this interpretation is no longer valid. Neverthe¬ 
less, in our analysis we used Eq. [2]as an empirical fitting 
function for data taken at all experimental pressures. In 
our case, the sum included all those bonds as in the ap¬ 
proach of Ref. The corresponding labeling scheme is 
shown in Eig. [3l 

The final ingredient of the model for 5'(Q,cj) is an 
empirical dispersion relation which was parameterized at 
all pressures as)i^ 

(TiuJcif = BqE Bh cos(27rh) + Bi cos(27r/) (3) 

+ Bhi[cos{27r{h -t- /)) + cos(27r(h — /))] 

+ B 2 h cos(47rh) + B 21 cos(47r/). 

Once again, we emphasize that our analysis was applied 
only to data in the (h,0,/) plane. Correspondingly, the 
above expression does not include any dispersion along 
b\ 

Eor a direct comparison with the measured neutron in¬ 
tensities, the model structure factor given in Eq. o was 
scaled with the magnetic form factor for Cu^+ written 
in the dipole approximation, and numerically convoluted 
with the resolution function of our measurement. The 
latter was a direct product of Gaussian resolution func¬ 
tions for momentum transfer in the principle h, I scat¬ 
tering plane and energy transfer, correspondingly. The 
instrumental wave vector resolution was determined by 
a Gaussian fit to the width of the nuclear Bragg reflec¬ 
tions in PHGG at very fine binning levels and was found 
to be practically isotropic (in the scattering plane) with 
(Jq = 0.054(2) A“^. The instrumental energy resolution 

= 0.069(3) meV was obtained by fitting the elastic 
line at a very fine binning level. 

Eor each pressure, the model was fitted simultaneously 
for the two two-dimensional data slices shown in Eig. [H 
At 9 and 18 kbar it was necessary to hold the parameter 
Bhi fixed at the value fitted at ambient pressure. At 9 
and 18 kbar, the upper part of the dispersion appears 
somewhat noisy and smeared out. The fits of the inten¬ 
sity (but not of the dispersion) therefore excluded the 
region above 2 meV to achieve reliable convergence. Ta¬ 
ble U shows the resulting empirical fit parameters for the 
dispersion relation as well as the resulting spin gap and 
band widths derived from these values. Overall, our re¬ 
sults at ambient pressure are in good agreement with the 
values obtained previously without the complication of 
performing an experiment inside of a pressure cell.— Our 
main result is that at 9 kbar, the gap A = 0.0(1) meV, 
i.e. the spectrum is gapless. Increasing the pressure to 
18 kbar, the gap remains closed with A = 0.0(1) meV. 


Table [III shows the fitted coefficients Ed at different 
pressures. Bonds 2-6 only exhibit modest changes in the 
whole 0-18 kbar region. The intensity coefficient for bond 
1 however, changes dramatically between 0 and 9 kbar 
from El = —1.11(5) meV at ambient pressure to Ei = 
—0.47(3) meV at 9 kbar but changes only moderately to 
El = —0.51(2) meV at 18 kbar. 


Parameter 

Stone et a/.— 

0 kbar 

9 kbar 

18 kbar 

Bo 

5.44(2) 

5.39(1) 

4.29(1) 

4.30(1) 

Bh 

2.06(3) 

1.76(1) 

2.0(2) 

2.0(2) 

Bi 

1.07(3) 

0.85(1) 

0.72(1) 

0.59(1) 

Bhi 

-0.39(1) 

-0.39(5) 

-0.39 

-0.39 

B 2 h 

-0.34(3) 

-0.55(2) 

-0.4(2) 

-0.5(1) 

B 21 

-0.22(2) 

-0.35(1) 

-0.36(1) 

-0.40(1) 

A 

0.98(6) 

1.05(3) 

0.0(1) 

0.0(1) 

Wh 

1.60(8) 

1.44(4) 

2.4(1) 

2.4(1) 

Wl 

1.18(8) 

1.04(4) 

1.7(1) 

1.7(1) 


TABLE I. Top: empirical dispersion parameters in meV^ at 
different pressures. Bottom: resulting spin gap A and band 
width Wh (wi) along h (1) in meV. 


Bond energy 

Stone et a/.— 

0 kbar 

9 kbar 

18 kbar 

El 

-1.3(3) 

-1.11(5) 

-0.47(3) 

-0.51(2) 

E 2 

0.7(3) 

0.87(6) 

0.77(4) 

0.75(3) 

E 3 

-0.3(1) 

-0.20(4) 

-0.14(3) 

-0.18(3) 

E 4 

0.1(3) 

0.23(5) 

0.22(3) 

0.29(3) 

Es 

-0.0(3) 

0.20(5) 

-0.08(3) 

-0.09(3) 

Ee 

-0.92(5) 

-0.85(4) 

-0.89(2) 

-0.89(2) 


TABLE 11. Fitted intensity modulation coefficients Ed in 
meV, for different pressures. The bond labeling corresponds 
to that in Fig. (3] 


V. DISCUSSION 

Our results demonstrate unambiguously that the mag¬ 
netic excitations in PHGG are gapless at and above 
9 kbar applied pressure. This is fully in agreement with 
the observation of long-range magnetic order above a 
QCP at Pc ~ 4.3 kbar using A pressure-induced 

reduction of the gap is also plausible in the context of the 
observation that Br-substitution, which leads to an in¬ 
crease of the lattice parameters and hence negative chem¬ 
ical pressure, causes the gap to increase/^ The discrep¬ 
ancy with Ref. ini where the gap seemed to extrapolate 
to zero only at much larger pressures, warrants an ex¬ 
planation. In that study, PHGG was assumed to possess 
a two-dimensional magnon dispersion. Under this as¬ 
sumption, to improve statistics, the data were integrated 
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along the /^-direction in a substantial range. As discussed 
above, the assumption is clearly invalid. Due to the dis¬ 
persion along 6*, any integration will produce an intensity 
maximum at the saddle-point energy (around 0.7 meV). 
We suggest that it is this saddle-point that was mistaken 
for the energy gap in Ref. U 

As stated above, for p > pc^ in the ordered state, the 
parameters E'd cannot be directly interpreted as bond en¬ 
ergies due to the anisotropic nature of the ordered state. 
This said, the fact that our isotropic model gives rea¬ 
sonably good fits to the experimental data in all cases, 
suggests that the anisotropy of spin correlations at high 
pressures remains small. This, in turn, suggests that the 
ordered moment remains small compared to the classical 
saturation value. In this case, a consistent and very pro¬ 
nounced decrease of \Ei\ with increasing pressure, and a 
lack of any drastic pressure dependence for the other pa¬ 
rameters in Table [3l are evidence that the transition is, 
in fact, driven by the weakening of a sin gle bond. That 
bond, labeled as “1” in Fig. [3] and in Rnf. llOl. happens to 
have the largest bond energy at ambient pressure. Fur¬ 
thermore, we note that we have not been able to observe a 
magnetic Bragg peak in the magnetically ordered phases 
which would be consistent with a small ordered moment 
size. 

We recall that there is some experimental /i+SR ev¬ 
idence that the magnetically-ordered state emerging at 
Pc is incommensurateOur wave vector resolution is 
approximately 0.07, 0.2 and 0.05 reciprocal lattice units 
FWHM along h, /c, and /, respectively. This sets an upper 
bound on the magnitude of the incommensurate propa¬ 
gation vector in the system. Beyond that, the present 
study is unable to confirm or disprove the incommensu¬ 


rate nature of the pressure-induced magnetic ordering. 
As far as the fitted dispersion and bond coefficients 
are concerned, only moderate changes are observed be¬ 
tween 9 kbar (in the first phase detected with /i+SR) and 
18 kbar (in the second phase found with //+SR). 


VI. CONCLUSION 

In summary, it appears that the magnetically ordered 
state induced in PHCC by hydrostatic pressure in excess 
of Pc ^ 4.3 kbar is gapless. The transition itself is driven 
primarily by the weakening of one particular exchange 
pathway. At high pressures PHCC is by no means a 
two-dimensional system, with a substantial magnon dis¬ 
persion along the third (6*) direction. 
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